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Abstract 

Gauge link smearing is widely used in lattice QCD computations. The 
idea is to remove the local (UV) fluctuations of the gauge field configu- 
rations while keeping the longer-range (IR) properties intact. Important 
applications are in the definitions of interpolating hadron operators as 
well as in updating the gauge field configurations with dynamical quarks 
in the so-called Hybrid Monte Carlo (HMC) algorithm. Here we study 
the effectiveness of various smearing methods and try to quantify these 
with tools known from Monte Carlo Renormalization Group (MCRG). 



1 Motivation 

In lattice gauge theory the group valued link variables encode the gauge field. 
In the Monte Carlo (MC) approach the space of gauge configurations is sampled 
according to the probability weight factor exp {—S), where S denotes the effec- 
tive gauge action, which may include terms coming from the integration of the 
fermionic variables as well. All observables are then determined from averaging 
functions of the gauge link variables over the MC-generated ensemble of gauge 
configurations. Wilson loops (traces of ordered products of gauge link variables 
along a loop) are the best known of such observables and allow to estimate 
the string tension. However, also quark propagators (leading to hadron n-point 
functions) are computed from the link variables. 

Since in many cases one is interested in the infrared behavior, i.e., correla- 
tions over distances large in lattice units, it is often worthwhile to smoothen the 
gauge configuration in a gauge-invariant fashion. This is done by constructing 
new link variables from an averaged sum of gauge link products along paths 
connecting the end points of the original link. Since the lattice size is not 
changed by this procedure we can consider such a smearing step as a MC block 
spin transformation (BST) with scaling factor 1. Actually, if we would be able 
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to determine the effective action corresponding to ttie smeared configurations 
weigfit factor, one could simulate directly the smoother gauge system, without 
changing long distance correlations. 

Smearing involves neighboring gauge links and thus has a certain effective 
radius of impact in lattice units. As long as the number of smearing steps 
is small compared to the lattice size, the long distance correlations will be 
unaffected. Iterating the smearing procedure extends this range and eventually 
the configurations will approach trivial ones, where all physical information is 
washed out. Non-local properties like topological modes will survive more of 
such steps and indeed smearing has also been used to improve identification of 
topological sectors. 

Here we will study the behavior of some observables under various smearing 
BSTs with the tools of the Monte Carlo Renormalization Group (MCRG) and 
try to obtain a more quantitative characterization of the smearing fiow and 
efficiency of these smearing methods. 

For completeness we summarize the four smearing algorithms that we have 
studied. 



APE Smearing [T]: The "smeared" link variable is built from a weighted 
sum over "staples" and the original link variable. 



O, X — ^ 



(1) 



Cf,u{n) = U^{n)Uf,{n + i))U},{n + fi) 

+ Ul{n-v)Uf,{n-v)U^{n-v + fi). (2) 



This sum is projected into X gSU(3) by maximizing Re tr[X[7^(n)t]. We take 
a = 0.55 as discussed in [2 . 

Hypercubic (HYP) Smearing [3j: This is a 3-step procedure which takes 
into account contributions from within a hypercube, again projected to SU(3). 
For the parameters involved in the definition we choose the values ai = 0.75, 
a2 = 0.6, a2 = 0.3, suggested in 

nHYP Smearing [4]: This is like HYP smearing, with a projection into U(3) 
via 

Vro^uidU']=U'^jl=,. (3) 

This smearing is differentiable with regard to the link variables, a property 
which is required for the HMC algorithm. The parameters are chosen like for 
the HYP smearing. 

STOUT Smearing ^5j: This also leads to a differentiable expression. The 
new link variable 

u'^{n)=eMiQ^{x))U^{n) (4) 

is constructed from the hermitian, traceless matrix which is built from a 
weighted sum over staples C^(n). This type of smearing is also widely used in 
simulations with the HMC algorithm. Following the notation of [S] we choose 
p = 0.1 as suggested in ^. 
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Figure 1: L.h.s.: The "step smearing function" for the plaquette observable 
Si, showing the mapping Si{ns) Si{ns + 1), where Hs denotes the number of 
smearing step (left- most point corresponds to =0). The points are connected 
by straight lines to guide the eye. R.h.s.: The subsector of 5*1 vs. 5*2 in the 
space of observables is shown to follow trajectories that are close to each other 
for the four smearing types considered. 

2 MCRG tools 

We study the different methods on an ensemble of 300 configurations of lattice 
size 16^ X 32, generated with two light dynamical quarks flavors and the Chirally 
Improved Dirac operator [71 [S] . The parameters of the action correspond to a 
lattice unit of 0.144 fm and a pion mass of 322 MeV. Details on action and 
simulation can be found in [H [TU] . 

The expectation values of five observables Si are determined: plaquette 5*1, 
planar rectangular 2 x 1-loop 5*2, bent loop S3, twisted bent loop ("chair") ^4 
and the topological charge S5 and susceptibility (in the formulation of [2 ). 

Figure [T] (l.h.s.) shows the change of the plaquette expectation value from 
one step of smearing to the next (for all four methods). Already here one 
finds a ratio of HYP : APE : STOUT w 3 : 2 : 1 (i.e., one HYP step corresponds 
approximately to three STOUT steps). The other observables show similar 
behavior. 

We note that for given configurations different smearing methods lead to 
different (close to integer) values of the topological charge for that configuration, 
stabilizing at a value typically only after 50-100 smearing steps. This nonlocal 
quantity thus proved not very useful for the subsequently discussed fiow study. 

In MCRG block spin transformations [11] the RG fiow in the space of cou- 
plings approaches a renormalized trajectory connecting fixed points. This renor- 
malized trajectory has a correspondence in the space of observables. Similarly, 
the smearing flow can also be followed in the space of observables (figure [U 
r.h.s.). These trajectories (also for other projections in the space of observ- 
ables) are close to each other for the different smearing types, while the speed 
along the trajectories can be different. We also find that HYP and uHYP are 
almost indistinguishable. 

In MCRG the block spin transformation of an ensemble of spin configura- 
tions distributed according to an action leads J2il^i^i to another ensemble with 
a different effective action '^i^'^Si, thus following a renormalized trajectory 
(connecting fixed points) in the space of coupling constants (/3i). We use this 
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Figure 2: The figures show the eigenvalues of ^("+1.") (nHYP results are close 
to HYP). 



framework to study the change of the effective action under smearing transfor- 
mations. The derivative matrix 



T' 



(n+l.n) 



d(3 



jk 



(5) 



(where the superscripts denote the number of blocking steps) defines a lineariza- 
tion of the flow in coupling space. It may be derived from cross-correlations of 
observables measured on conflgurations for different numbers of smearing steps 
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For BSTs close to a fixed point approaches the fixed point value T* and 

then its eigenvalues Xi (if greater than 1) can be related to the relevant critical 
exponents through I'i — ln&/lnAi (here b is the scale factor of the EST). 

In our case we are interested in ^("=+1."=) for small Ug- Its diagonalization 
gives the information about the main fiow directions and the speed of fiow. 
Actually, in case one knows the starting point in coupling space one can 
approximately recover [Pi)' . Related problems have been discussed recently in 
the context of trivializing maps [13] . For the gauge configuration ensembles used 
here we do not know (/3i), since they result from a simulation with dynamical 
fermions and thus the full effective gauge action is not the simple one used for 
the gauge part alone. 

Figure [2] shows the leading eigenvalues of the diagonalization. Only two 
of them are significantly larger than 1; asymptotically they all approach 1 of 
course. The data for HYP smearing show a peak in the second step, they also 
give the largest eigenvalues indicating the strongest effect (on the observables 
considered). 

Comparing the eigenvalues with the (normalized) eigenvectors for the leading 
eigenvalue (figure [3]) we find that the dominant direction is the plaquette Si 
followed by the bent loop in all cases. The eigenvectors are defined only up to 
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Figure 3: Here the eigenvector components of the leading eigenvalue of T*^""*"^'") 
are plotted (nHYP results are not shown here but are close to HYP). 



normalization (we normalized to unit length) and overall sign, which we choose 
such that continuity in Us is obtained. 

The eigenvectors for the 2nd largest eigenvalue confirms that APE and 
STOUT smearing follow similar trajectories. Here the bent loop and the twisted 
bent loop dominate the fiow. The eigenvectors for HYP smearing exhibit 
stronger fluctuation: more directions in coupling space are activated, mixing 
the local degrees of freedom stronger than the other types. 

3 Conclusions 

We try to make a step towards a systematic classification of gauge link smearing 
algorithms. Viewing smearing steps as analogues of block spin transformations, 
we attempt to quantify their effectiveness. This is a first attempt towards devel- 
oping a tool for classification of the effect smearing has in configuration space. 
We find that the different methods follow similar trajectories in the space of 
local observables and effective couplings, however with different speed. The 
nHYP results are close to the HYP results. In our study only few local observ- 
ables have been included; considering more would allow for a better and more 
complete quantification. 
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